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ABSTRACT 

We present photometric and spectroscopic observations of the interacting transient SN 2009ip 
taken during the 2013 and 2014 observing seasons. We characterise the photometric evolution 
as a steady and smooth decline in all bands, with a decline rate that is slower than expected 
for a solely 56 Co-powered supernova at late phases. No further outbursts or eruptions were 
seen over a two year period from 2012 December until 2014 December. SN 2009ip remains 
brighter than its historic minimum from pre-discovery images. Spectroscopically, SN 2009ip 
continues to be dominated by strong, narrow (<2000 kms -1 ) emission lines of H, He, Ca, 
and Fe. While we make tenuous detections of [Fe n] 47155 and [O i] 446300,6364 lines 
at the end of 2013 June and the start of 2013 October respectively, we see no strong broad 
nebular emission lines that could point to a core-collapse origin. In general, the lines appear 
relatively symmetric, with the exception of our final spectrum in 2014 May, when we observe 
the appearance of a redshifted shoulder of emission at +550 k m s'. The lines are not blue- 
shifted, and we see no significant near- or mid-infrared excess. From the spectroscopic and 
photometric evolution of SN 2009ip until 820 days after the start of the 2012a event, we 
still see no conclusive evidence for core-collapse, although whether any such signs could be 
masked by ongoing interaction is unclear. 

Key words: stars: massive — supernovae: general — stars: mass loss - supernovae: individ¬ 
ual (SN2009ip) 


1 INTRODUCTION 

The physical distinction between Type Iln supernovae (SNe) and 
SN impostors is clear: Type Iln SNe arise from massive star£] 
undergoing core-collapse within a dense circumstellar medium 
(CSM), while SN impostors are eruptions or outbursts from mas¬ 
sive stars which do not undergo core-collapse l lVan Dvk et al.l200(il . 
and references therein). Making this distinction from an observa¬ 
tional basis is more difficult, as for both Type Iln SNe and SN im¬ 
postors, the evolution and properties are influenced by the interac¬ 
tion of fast ejecta with the CSM, with similar charac teristic nar- 
row line emission in the observed spectra (e.g. lMaund et al.ll2006l : 
ISmith et al .11201 1 1) . 

One of the best studied and most controve rsial SN impostor s 
to date is SN 2009ip. First identified in 2009 bv lMaza et alJd2009t ), 
and initially given a supernova designatior^}, it was re-classified 
as an impostor on the basis of its narrow emission line s, low ab¬ 
solute magnitude (V — 14) and previous variability (Ber ger et al.l 
1 20091 : iNTiller et al.l2009l) . Pre-explosion images of NGC 7259 taken 
with the Hubble Space Telescope ( HST) revealed the progenitor of 
SN 2009ip to be a luminous ( V —10. 3) star with an estimated zero- 
age m ain sequence mass of ~60 M 0 dSmith et al.ll201(jl ; lFolev et al.l 
1201 lh . SN 2009ip exhibited sporadic variability over 3 years after 
discovery, with multiple outbur sts reaching absolut e magnitudes 
between R =-12 and -14 mag dPastorello et al.1120131) . During this 
period, the spectra were dominated by ~700-1300 km s -1 em ission 
lines of H, along with He, Ca, Na a nd Fe dPastorello et al.l[20131 : 
ISmith et all 1 201 ()l : iFolev et alj|201lh . Intriguingly, there was also 
absorptio n at v e locities o f up to 12500 kms -1 as early as Septem¬ 
ber 2011 dPastorello et alj|2013l) . implying that at least some ma¬ 
terial had been ejected at velocities more usually associated with 
core-collapse SNe. 


* E-mail:mf@ast.cam.ac.uk 

1 While some interacting transients have been linked to thermonuclear ex¬ 
plosions iHamuv et al.l2003h . we will not discuss these here. 

2 In this paper, we use the designation ‘'SN 2009ip” to refer to the object 
in question’s outburst in 2009, the subsequent period of variability, and the 
2012a and 2012b events. 


In 2012 August, SN 2009ip underwent a new eruption lasting 


dPrieto etal. 2013; Pastorello et al 

2013); Mauerhan et al. 2013}; 

Graham et al.. 20141; Margutti et al. 

2014). This eruption was fol- 


lowed by a second outburst ~40 days later, which peaked at an 
absolute magnitude of R =-18 mag. These two re-brightenings of 
SN 2009ip were termed the “2012a” and “2012b” events respec¬ 
tively. The 2012a event showed a ~20 day rise in all optical bands, 
followed by a decline on a similar timescale. Broad (~8000 km s -1 
minima) P Cygni lines were seen in H, He and Na, together with 
strong, narrow emission lines above a blue continuum. Immedi¬ 
ately after the 2012a event, SN 2009ip re-brightened dramatically 
during the 2012b event, increasing by nearly 4 magnitudes in two 
weeks dPrieto et~ahll 20131) . The broad absorptions in the spectra of 
SN 2009ip disappeared at the start of the 2012b event, leaving a 
blue continuum with narrow, Lorentzian profile lines. Following 
the maximum of the 2012b event, SN 2009ip underwent a slower 
decline, which was accompanied by the reappearance of broad ab¬ 
sorptions witfficomffiexjine profiles and multiple absorption com¬ 
ponen ts |Fraser et al.ll2013al; Pastorello et alJl2oi~3l : iMareutti et al.l 
|20I4|) . iFraser et al] and lMargutti et al.l usedthe magnitude of SN 
2009ip at the end of the 2012b event to place a limit of <0.02 M 0 
on the mass of any radioactive 56 Ni synthesised by SN 2009ip. 


The CSM surrounding SN 2009ip clearly plays a significant 
role in the evolution of t he transient, and h a s been studied in de¬ 
tail by multipl e groups 1 Ofek et al] l2013bl : iMargutti et alj 120141 ; 
IFox et al.l2Q15l) . lofek et al.1 J2013bl) reviewed all the available mass 
loss diagnostics for SN 2009ip, including X-rays, radio, and the 
narrow Balmer lines, and concluded that the average mass-loss rate 
for the progenitor over the years prior to explosion was between 
10~ 2 and 10~ 3 M 0 per year. This value was lower tha n the mass 
loss inferred from the Her luminosity of SN 2009ip dOfek et afl 
l2013bl : IFraser et akl l2013al) . but according to the authors, could 
be reconciled with it by assuming a non-spherical distribution of 
CSM. Fro m a d etailed study of the Balmer decrement in SN 2009ip, 
iLevesque et al.l d20l4 ) argued that the CSM was likely aspherical, 
and possibly involved a disk-like structure; wh ile spe ctropolarimet- 
ric observations presented bv lMauerhan et all d2014l) also indicated 
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a complex geometry for SN 2009ip. 'Mauerhan et al.1 suggested that 
the narrow emission for SN 2009ip arose from an inclined disk, 
while the broad components came from a distinct, and perpendicu¬ 
lar emitting region. 

After the 2012a and 2012b events, SN 2009ip appears to be 
continuing to fade in all bands. The decline rate of SN 2009ip 
slowed from that initially observed after the peak of the 2012b 
event, and as of mid 2013 appeared to be approximately half that 
expected from the decay of radioactive 56 Co cGraham et alj|2014h . 


There have been no further reported outbursts or re-brightenings 
(Grah am et alj) . however, the magnitude of SN 200 9ip still re¬ 
main s brighter than tha t of the presumed progenitor I Sm ith et~aD 
|20 1(1 [Folev et al.ll201ll) . The late time (post-2012b) spectra of SN 
2009ip reveal relatively n arrow (~l-2x!0 3 km s -1 ) emission from 
H, Fie, Na, Ca, and Fe dGraham et ai] |20l4 ISmithetalJl20Tl 
Margutti et al. 2014), albeit without any of the strong nebular lines 
typically associated with core-collapse SNe. We note, however, that 
Ifox et all I2OI5I) found evidence for very weak [O 1 ] emission in a 
spectrum obtained at late phases. 

Despite the rich observational dataset obtained for S N 2009ip 
acros s an exceptionally wide range of wavelengths l lMargutti et al.1 
i2014h . the interpretation of the events of late 2012 remains uncer¬ 
tain. During the 2012a and 2012b events, ejecta were seen at ve¬ 
locities of over 10 4 kms -1 which together with the unprecedented 
luminosity, was suggested by iMauerhan et al.l d20l 3l) to be evi¬ 
dence that SN 2009ip had undergone core-collapse and was now 
a genuine Type Iln SN. [Smith et aTI |2014|) expanded on this sce¬ 
nario, and proposed that the 2012a event was the core-collapse of 
a blue supergiant, while the 2012b event was powered by the sub¬ 
sequent collision of t he ejecta from this event with the surround- 
ing CSM. In contras t. IPastorello et alJ d2013l) . lFraser et al] d2013al) 
and iMargutti et alj 1 2014) remained open as to the final fate of 
SN 2009ip. Pastorello et al.1 argued that as high velocity absorp¬ 
tion had already been observed in 2011, this could not be taken as 
proof positive of core-collapse during either of the 2012a or 2012b 
events. The absence of emission lines associated with nucleosyn- 
thesised elements, and the relatively restrictive limit to the ejected 
56 Ni mass could also point to a non-terminal eruption. If not pow¬ 
ered by a core-collapse. iFraser et al.l and [Margutti et al.1 suggested 
that the luminosity of the 2012b event might be explained by the 
efficient conversion of the kinetic energy of a relatively small erup¬ 
tion as it collided with t he CSM, An alternativ e explanation for SN 
2009i p was explored bvISoker & KashU d2013l) . lKashi et all d2013l) 
and lTsebrenko & Sokeil ( 2013 ), who proposed that SN 2009ip may 
be a “mergerburst” resulting from two massive stars in a binary sys¬ 
tem. In such a scenario, the secondary may either merge completely 
with the primary, or accrete a sign ificant fraction of its mass onto 
the primary during periastron dKashi et alj2013ll . The high velocity 
material seen in SN 2009ip could then be formed by the interaction 
of a fast jet of material, launched durin g the merger, with a cir- 
cumbinary envelope dTsebrenko & SokenuOl 3I) . 

Regardless of whether the 2012a or 2012b events of SN 
2009ip were terminal or non-terminal events, the mechanism be¬ 
hind the outbursts seen prior to and including late 2012 remain un¬ 
known. There have been several theoretical proposals for increased 
mass loss shortly before core-collapse c aused by vigorous convec¬ 
tion during C and la ter burning stages JOuataert & Shiodal2012l : 
ISmith & Arnettl20l4) . Another candidate m echanism is the “pulsa- 
tional pair instability” dWooslev et alJl2007h . which can lead to the 
release of a large amount of energy, and the ejection of several so¬ 
lar masses of materia]_without destroying the progenitor star. Alter¬ 
natively, [Soker & Kashil d2013l) proposed that SN 2009ip could be 


explained by a mergerburst, where a ~60 M 0 primary and a lower 
mass main-sequence secondary merged. In this scenario, the previ¬ 
ous outbursts of SN 2009ip during the period 2009 to 2012 August 
may have been triggered by close periastron passages before the 
final merger event. 

The Public ESO Spectroscopic Survey for Transient Objects 
(PESSTO: ISmartt et alj20l3) observation campaign for SN 2009ip 
concluded in late 2012 De cember. In 2013 Ap ril, SN 2009ip was 
visible again from La Silla ( IFraser et al.l2013g) : this paper presents 
observations of SN 2009ip taken during the 2013 observing sea¬ 
son until 2013 December, along with observations taken between 
2014 April and December. In particular, we address the questions 
of whether there is new evidence for core-collapse, if there are sig¬ 
natures of dust formation in the ejecta or shocked CSM gas, and 
crucially, whether SN 2009ip has survived the events of late 2012. 
The true fate of SN 2009ip is of considerable interest, especially 
given the rec ent identification of pre-explosion outbursts in sev¬ 
eral SNe (e.glPastorello et alJl2007l : lOfek et ai1l2013al : IFraser et all 


12013d : lOfek et all 2014h . SN 2009ip is also one of the best moni¬ 
tored SNe of any type, with 5 years of data since its discovery, and 
as such represents a unique window into the lives of massive stars. 

In all of the following we adopt a distance modulus of 
31.55 mag, and a foreground extinction of A«f=0.051 mag to¬ 
wards SN 2009i p, con sistent with that used bv ISmith et al.l d2010h : 
IPastorello et alj d2013h : IFraser et alj d2013ah . We also adopt 2012 
September 23 (MJD 56193) as the start of the 2012b event, all 
phases quoted in the paper are with respect to this date. 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Imaging 

Optical imaging of SN 2009ip was obtained with the 3.6m New 
Technology Telescope (NTT) + ESO Faint Object Spectrograph 
and Camera 2 (EFOSC2), the 1.3 m Small and Moderate Aperture 
Research Telescope System (SMARTS) telescope + ANDICAM 
(operated by the SMARTS Consortium), the 2 m Faulkes South 
telescope + Spectral imager (FS03) and the 2 m Liverpool Tele¬ 
scope + Infrared-OpticaLOptical (10:0) camera. The EFOSC2 im¬ 
ages were reduced using the PESSTO pipeline, as described in 
ISmartt et al.l l20l4 . while the ANDICAM, FS03 and 10:0 im¬ 
ages were automatically reduced with their associated instrument 
pipelines. In addition to these, we used the filtered acquisition im¬ 
ages taken before our Very Large Telescope (VLT) + XShooter and 
Gemini South + Gemini Multi-Object Spectrograph (GMOS) spec¬ 
tra. The XShooter acquisition image is taken with a separate acqui¬ 
sition camera to the main instrument, for which biases and flat- 
fields are not available. However, the image did not show a strong 
gradient in the background, and so we assume the flat field correc¬ 
tion is small. The Gemini images were reduced using the tasks in 
the gemini package within iraf. The ANDICAM, 10:0, FS03 and 
EFOSC2 images were processed with lacosmic to remove cosmic 
rays, before consecutive exposures taken with the same instrument 
and filter were coadded. As only single, short exposures were taken 
with GMOS and the XShooter acquisition camera, these data were 
not coadded or processed with lacosmic. 

For all images, point-spread function (PSF) fitting photome¬ 
try was carried out using snoopy, a custom iraf package based on 
daophot. The PSF was modelled using point sources on each frame, 
and was fit simultaneously to both SN 2009ip and the nearby fore¬ 
ground star at (~5"to the NE at 22 h 23 m 5!51 -28°56'4779) to ac¬ 
count for any additional flux from the latter. The fit was iterated 
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to minimise the residual, while artificial star tests at the SN posi¬ 
tion were used to estimate the uncertainty in magnitude. The ze- 
ropoint for each frame was determined u sing aperture ph otome- 
try of the tertiary sequence stars listed by IPastorello et all l2013t ), 
with the uncertainty in the zeropoint taken to be the standard devi¬ 
ation of the values calculated from each of the individual sequence 
stars. Aperture corrections were also measured from the sequence 
stars, and where necessary applied to the photometry of SN 2009ip, 
and c olou r terms were applied to the measured magnitudes as in 
iFraser et~ahl l2013ah . 

At these late epochs, the spectrum of SN 2009ip is dominated 
by line rather than continuum emission. As a consequence of this, 
the photometry of SN 2009ip is sensitive to the precise transmis¬ 
sion of each filter used. As shown in Fig. [I] the Bessell and Sloan 
filters from the various instruments used for observations of SN 
2009ip have somewhat different transmissions. These differences 
result in flux from various lines being included or excluded from 
a particular passband - for example the EFOSC2 i filter includes 
flux from the T8662 line of the Ca II NIR triplet, while the 10:0 i 
filter does not. As a check on the significance of the flux present in 
line emission at this epoch, we took the CMOS spectrum from Au¬ 
gust 12 and set the flux in the region of Ha- to the continuum level. 
Synthetic Bessell R-band photometry performed with iraf synphot 
on both the original and Hu-free spectra gives a magnitude for the 
latter which is 1.3 mag fainter, demonstrating line flux can have a 
significant effect on our measured photometry. 

It is possible to ameliorate this to some extent using_synthetic 
photometry of spectra (“S-correction”, e.g. lStritzinger et alil2002i) , 
however this is dependent on the availability of carefully measured 
total transmission functions (i.e. telescope optics, filter and CCD), 
which were unavailable for several of the instruments we have used. 
Hence we have not performed a full S-correction for all filters, but 
rather have applied a simple correction using the available filter 
functions, for the instruments and filters where there was the great¬ 
est disagreement. We have also taken the EFO SC2 fi lters as our 
reference filter system, rather than the standard feesse 11 d 199(1) filter 
curves, with the motivation of both reducing the number of correc¬ 
tions which we make to our data (i.e. the EF0SC2 data does not 
have to be corrected), and also to ensure consistency with the data 
presented in lFraser et al.1 d2013ah . Average colour ter ms to convert 
EFOS C2 photometry to other systems are available in Sm artt et~ahl 
d2014l) . 

For each photometric point to be corrected, we have taken the 
closest available spectrum of SN 2009ip, and with synphot, mea¬ 
sured the difference in synthetic magnitude when using the specific 
filter function for that instrument, and when using the correspond¬ 
ing EF0SC2 filter (AAf^_ £f0SC2 ). This difference cannot be ap¬ 
plied directly to the SN 2009ip magnitudes, as some fraction of it 
will already be accounted for by the applied colour term. To esti¬ 
mate the portion of this difference which is accounted for by the 
colour term, we first estimated the effective temperature of the SN, 
using a blackbody fit to the spectrum (excluding the regions around 
Ho- and H/f). We then took a Kurucz model at a similar tempera¬ 
ture, and appropriate to a solar metallicity dwarf, on which we also 
measured the difference in synthetic magnitude (AM 2MS _£ f0SC2 ) 
using the filter functions for the instrument in question, and the 
filter function for EFOSC2. This difference is what one would ex¬ 
pect for a continuum-dominated source, such as those stars from 
which the applied colour terms are derived. Finally, we can cal¬ 
culate (A Mf»s -EFOsa) {^^ins-efoscz)' w hich is the additional 
correction to the measured photometry of the SN at that epoch (in 
addition to the colour term). The final calibrated magnitudes for SN 



4000 5000 6000 7000 8000 9000 
Observer frame wavelength (A) 


Figure 1. Filter functions for the filters used for the followup of SN 2009ip; 
for FS03 these are transmission functions which also include the contri¬ 
bution of the Earths atmosphere. All transmission filters have been nor¬ 
malised to 100% transmission at their peak. For comparison, an observer 
frame spectrum of SN 2009ip from 2013 April is shown in grey. 


2009ip in the EFOSC2 (~Landolt) system, together with the total 
uncertainty (comprised of the error in SN magnitude and zeropoint, 
added in quadrature) are reported in TableQ] For the LT+IO:0 im¬ 
ages which were taken in Sloan-like r and /-filters, we were unable 
to homogenise the measured photometry with the data from other 
telescopes, even after applying the correction from synthetic pho¬ 
tometry. We have hence excluded these images from our analysis. 
With the exception of the LT ri images, the measured magnitudes 
for SN 2009ip using different telescopes on the same night are con¬ 
sistent to within the uncertainties. 

We also took near-infrared (NIR) imaging of SN 2009ip in 
JHKs filters using both the NTT + Son of ISAAC (SOFI) and the 
2.6m Nordic Optical Telescope + NOTCam. The SOFI data were 
reduced using the PESSTO pipeline dSmartt et alj|2014l) , while the 
NOTCam images were reduced using the external notcam pack¬ 
age (version 2.5) within iraf. In both cases, the images were flat- 
fielded, and sky-subtracted using a sky frame constructed from 
median combined on-source dithers. PSF-fitting photometry was 
performed on the SN using the same procedure as for the optical 
imaging, although the zeropoint was determined directly from cat¬ 
alogued Two Micron All-Sky Survey (2MASS) sources in the field, 
and no colour terms were applied. The calibrated magnitudes are 
reported in Table[2] For epochs where SN 2009ip could not be seen 
in these images, a limiting magnitude was estimated from artificial 
star tests at the expected position of the transient. 

Mid-infrared (MIR) imaging was taken of SN 2009ip at mul¬ 
tiple epochs in 2013 with the Spitzer Space Telescope + Infrared 
Array Camera (IRAC), as listed in Table [2] As Spitzer is now in 
its warm mission phase, SN 2009ip was observed only in Chan¬ 
nel 1 (3.6qm) and Channel 2 (4.5qm). Images from the Spitzer 
archive taken in 2009 and 2010 were also examined, but the posi¬ 
tion of the transient lies off the edge of the detector in these data, 
and hence these were not considered any further. Photometry was 
performed on the post-Basic Calibrated Data (pBCD) images us- 
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Table 1. Calibrated optical magnitudes of SN 2009ip in the Landolt system, with associated uncertainties in parentheses. Phase is with respect to the start of 
the 2012a event (MJD 56193). 


Phase (d) 

Date 

MJD 

B (mag) 

V (mag) 

R (mag) 

I (mag) 

Instrument 

192 

2013-04-03.40 

56385.40 

- 

19.23 (0.06) 

18.24 (0.06) 

- 

EFOSC2 

193 

2013-04-04.38 

56386.38 

- 

19.25 (0.05) 

- 

- 

EFOSC2 

195 

2013-04-06.41 

56388.41 

19.95 (0.10) 

- 

- 

- 

EFOSC2 

197 

2013-04-08.40 

56390.40 

- 

19.26 (0.10) 

18.32 (0.07) 

18.86 (0.08) 

ANDICAM 

201 

2013-04-12.39 

56394.39 

- 

19.43 (0.07) 

18.42 (0.08) 

18.76 (0.10) 

ANDICAM 

202 

2013-04-13.37 

56395.37 

- 

19.25 (0.06) 

- 

- 

EFOSC2 

203 

2013-04-14.37 

56396.37 

- 

19.34 (0.07) 

- 

- 

EFOSC2 

205 

2013-04-15.90 

56397.90 

- 

19.34 (0.11) 

18.29 (0.06) 

18.80 (0.09) 

ANDICAM 

209 

2013-04-20.37 

56402.37 

- 

19.38 (0.05) 

- 

- 

EFOSC2 

219 

2013-04-29.85 

56411.85 

- 

- 

18.47 (0.05) 

- 

ANDICAM 

222 

2013-05-03.37 

56415.37 

- 

- 

18.46 (0.12) 

- 

GMOS 

225 

2013-05-06.33 

56418.33 

- 

- 

18.42 (0.06) 

- 

XShooter 

239 

2013-05-20.34 

56432.34 

- 

- 

18.55 (0.06) 

- 

GMOS 

242 

2013-05-23.40 

56435.40 

- 

- 

18.59 (0.05) 

- 

GMOS 

245 

2013-05-25.78 

56437.78 

20.48 (0.09) 

19.89 (0.06) 

18.53 (0.05) 

18.96 (0.08) 

FS03 

248 

2013-05-28.79 

56440.79 

- 

19.93 (0.11) 

18.61 (0.07) 

- 

FS03 

254 

2013-06-03.69 

56446.69 

20.55 (0.05) 

20.02 (0.07) 

18.62 (0.04) 

19.12(0.08) 

FS03 

254 

2013-06-04.34 

56447.34 

- 

- 

18.67 (0.03) 

- 

GMOS 

264 

2013-06-14.35 

56457.35 

- 

- 

18.77 (0.05) 

- 

GMOS 

269 

2013-06-18.82 

56461.82 

20.34 (0.05) 

20.08 (0.05) 

18.70 (0.04) 

19.22 (0.04) 

FS03 

280 

2013-06-30.40 

56473.40 

- 

- 

18.82 (0.03) 

- 

GMOS 

293 

2013-07-12.81 

56485.81 

20.73 (0.04) 

20.17 (0.07) 

18.76 (0.04) 

19.34 (0.08) 

FS03 

303 

2013-07-23.33 

56496.33 

- 

- 

18.85 (0.05) 

- 

GMOS 

306 

2013-07-26.11 

56499.11 

20.83 (0.23) 

20.02 (0.10) 

- 

- 

10:0 

308 

2013-07-27.70 

56500.70 

20.67 (0.07) 

20.12(0.06) 

18.82 (0.04) 

- 

FS03 

308 

2013-07-28.34 

56501.34 

20.68 (0.08) 

20.06 (0.05) 

18.97 (0.06) 

19.27 (0.05) 

EFOSC2 

311 

2013-07-30.70 

56503.70 

20.91 (0.11) 

20.36 (0.10) 

18.89 (0.06) 

19.58(0.13) 

FS03 

316 

2013-08-05.20 

56509.20 

20.71 (0.06) 

20.08 (0.05) 

19.01 (0.06) 

19.45 (0.06) 

EFOSC2 

317 

2013-08-05.79 

56509.79 

20.70 (0.05) 

20.14(0.08) 

18.82 (0.04) 

19.44(0.10) 

FS03 

319 

2013-08-08.02 

56512.02 

20.67 (0.09) 

20.02 (0.09) 

- 

- 

10:0 

323 

2013-08-11.62 

56515.62 

20.78 (0.08) 

20.27 (0.07) 

18.90 (0.04) 

- 

FS03 

323 

2013-08-12.15 

56516.15 

- 

- 

18.93 (0.12) 

- 

GMOS 

326 

2013-08-15.13 

56519.13 

20.82 (0.04) 

20.20 (0.05) 

19.03 (0.07) 

19.52 (0.05) 

EFOSC2 

327 

2013-08-16.07 

56520.07 

20.88 (0.12) 

20.06 (0.05) 

- 

- 

10:0 

329 

2013-08-17.59 

56521.59 

20.88 (0.09) 

20.39 (0.05) 

18.96 (0.03) 

19.59 (0.10) 

FS03 

338 

2013-08-26.63 

56530.63 

- 

20.20 (0.19) 

18.99 (0.12) 

- 

FS03 

338 

2013-08-27.10 

56531.10 

20.88 (0.06) 

20.27 (0.06) 

19.14(0.06) 

19.61 (0.09) 

EFOSC2 

341 

2013-08-29.62 

56533.62 

- 

- 

19.01 (0.08) 

- 

FS03 

341 

2013-08-30.01 

56534.01 

20.96 (0.12) 

20.09 (0.07) 

- 

- 

10:0 

351 

2013-09-09.04 

56544.04 

- 

20.17 (0.07) 

- 

- 

10:0 

354 

2013-09-12.20 

56547.20 

20.75 (0.05) 

20.14(0.04) 

19.12(0.05) 

19.66 (0.08) 

EFOSC2 

361 

2013-09-18.95 

56553.95 

- 

20.24 (0.23) 

- 

- 

10:0 

369 

2013-09-26.95 

56561.95 

20.77 (0.08) 

20.26 (0.07) 

- 

- 

10:0 

376 

2013-10-04.05 

56569.05 

20.99 (0.06) 

20.39 (0.10) 

19.23 (0.04) 

19.77 (0.07) 

EFOSC2 

388 

2013-10-16.18 

56581.18 

20.86 (0.13) 

20.46(0.15) 

19.30 (0.06) 

19.82 (0.08) 

EFOSC2 

397 

2013-10-25.14 

56590.14 

21.07 (0.10) 

20.52 (0.11) 

19.34 (0.08) 

19.96 (0.24) 

EFOSC2 

409 

2013-11-05.83 

56601.83 

- 

20.41 (0.32) 

- 

- 

10:0 

427 

2013-11-23.87 

56619.87 

21.24 (0.09) 

20.51 (0.07) 

- 

- 

10:0 

430 

2013-11-26.82 

56622.82 

21.25 (0.13) 

20.50 (0.07) 

- 

- 

10:0 

434 

2013-12-01.05 

56627.05 

21.19(0.06) 

20.67 (0.05) 

19.43 (0.04) 

19.98 (0.09) 

EFOSC2 

457 

2013-12-24.05 

56650.05 

21.25 (0.07) 

20.71 (0.04) 

19.45 (0.05) 

20.13 (0.05) 

EFOSC2 

575 

2014-04-20.32 

56768.32 

- 

21.05 (0.05) 

19.80 (0.07) 

20.40 (0.27) 

EFOSC2 

577 

2014-04-22.36 

56770.36 

21.62 (0.10) 

- 

- 

- 

EFOSC2 

578 

2014-04-23.40 

56771.40 

- 

21.17 (0.06) 

19.84 (0.05) 

20.65 (0.10) 

EFOSC2 

583 

2014-04-28.38 

56776.38 

- 

21.03 (0.18) 

19.80 (0.11) 

20.52 (0.21) 

EFOSC2 

589 

2014-05-04.37 

56782.37 

21.61 (0.61) 

21.10(0.09) 

19.90 (0.08) 

20.83 (0.29) 

EFOSC2 

692 

2014-08-16.12 

56885.12 

21.92 (0.08) 

21.45 (0.05) 

20.21 (0.04) 

20.89 (0.08) 

EFOSC2 

723 

2014-09-16.20 

56916.20 

- 

- 

19.90 (0.08) 

- 

EFOSC2 

731 

2014-09-24.02 

56924.02 

21.79 (0.13) 

21.45 (0.05) 

20.27 (0.06) 

- 

EFOSC2 

736 

2014-09-29.21 

56929.21 

21.94 (0.10) 

21.54 (0.07) 

20.28 (0.05) 

- 

EFOSC2 

782 

2014-11-14.06 

56975.06 

21.99 (0.09) 

21.61 (0.09) 

20.44 (0.06) 

21.21 (0.08) 

EFOSC2 

821 

2014-12-23.06 

57014.06 

22.00 (0.07) 

21.73 (0.07) 

20.55 (0.07) 

21.45 (0.31) 

EFOSC2 
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Figure 2. MIR images of SN 2009ip, taken with the Spitzer Space Telescope 
+ IRAC on 2013 July 25. 


analog was observed at comparable airmass to SN 2009ip and was 
used to correct for the regions of low atmospheric transmission in 
the NIR. In the final XShooter spectrum obtained of SN 2009ip in 
2014 August, no flux from SN 2009ip was detected in the NIR arm. 

NIR spectroscopy was obtained with NTT+SOFI. The SOFI 
spectrum was taken with the Blue Grism (BG) and reduced using 
the PESSTO pipeline; the steps involved were flat-fielding, sky- 
subtraction using pairs of dithered observations and wavelength 
calibration. Water vapour absorption was identified and corrected 
for in the spectrum using a telluric standard (a solar analog) which 
was observed immediately after the spectrum of SN 2009ip, and 
at a similar airmass. An approximate flux calibration of the spec¬ 
trum was also made using the telluric standard. We discuss the NIR 
spectra obtained for SN 2009ip in Section l4~2l 


ing the mopex package llMakovoz & M arleau 20(3). A mosaic of 
the field was created (an example of which is shown in Fig. H, 
and sources were identified on this image, before Point Response 
Function (PRF) fitting photometry was performed on the individual 
frames. Measured fluxes were then corrected for normalisation as 
per the IRAC Instrument Handbook. No array-position dependent 
colour correction was made to the data, although as the SN is not 
particularly blue at these late epochs, the effect of this should be 
negligible. The measured fluxes for SN 2009ip are listed in Table 

H 

2.2 Spectroscopy 

We obtained a series of intermediate resolution optical spectra of 
SN 2009ip between 2013 April and December, together with a 
subsequent set of spectra from the 2014 observing season cov¬ 
ering the period 2014 April to September, as listed in Table [3] 
Longslit spectra were taken with the NTT+EFOSC2 (as part of the 
PESSTO project), the 8.2m Very Large Telescope (VLT) + FO- 
cal Reducer and Spectrograph 2 (FORS2), the 8. lm Gemini South 
telescope + GMOS2, and the 10m Keck II tel escop e + the DEep 
Imaging Multi-Object Spectrograph (DEIMOSi lFaber et alJi2003l ). 
Higher resolution echelle spectra were taken with VLT+ XShooter 
dVernet et alj201lh ; the XShooter spectra cover a broad wavelength 
range from the atmospheric cut off in the blue, to 2.5 pm. We also 
obtained three spectra with the ANU 2.3m telescope + Wide-Field 
Spectrograph (WiFeS), these had low S/N even after rebinning to a 
lower resolution, but are shown for completeness in Fig. [3] 

All EFO SC2 s pectra were reduced using the PESSTO pipeline 
JSmartt et alJl2014i) . while the FORS2 and GMOS spectra were re¬ 
duced within iraf using the longslit and gemini gmos packages re¬ 
spectively. In all cases, the data were reduced in the standard fash¬ 
ion: extracted, wavelength and flux calibrated, and in the case of the 
EFOSC2 spectra corrected for telluric absorption. The WiFeS spec¬ 
tra were reduced using the pywifes pipeline I Chi ldress et all2014h ; 
as WiFeS is an integral field unit spectrograph, one dimensional 
spectra were extracted from the data cube using a PSF-weighted 
fit. The XShooter spectra were reduced with the ESOREX recipes 
run under the reflex environment. For all XShooter spectra, a solar 

3 PESSTO data taken up to 2014 April 30 are available from the ESO 
archive as Phase III products (within the Spectroscopic Survey Data Re¬ 
lease 1; SSDR1), and the remainder of the data will be released in 2015 as 
SSDR2. Detailed instructions on obtaining PESSTO data is available on the 
PESSTO website, www.pessto.org. 


3 PHOTOMETRIC EVOLUTION 

The late time lightcurve of SN 2009ip is shown in Fig.[5] SN 2009ip 
appears to be slowly declining in magnitude, with no further out¬ 
bursts or re-brightenings over the two years following the 2012a 
and 2012b events. Extrapolating the lightcurve back to 2012 De¬ 
cember, we find that the slow decline is consistent with the last 
photometry taken at the end of 2012 December. Most notably, the 
decline rate is significantly slower than the 0.98 mag per 100 days 
which is typical of the late-time 56 Co-powered tail of a SN. From 
a linear fit to the R-band lightcurve in the period +200-450 days 
(relative to the start of the 2012a event), we measure a decline rate 
of 0.34 mag per 100 days, while the V-band yields a decline rate 
of 0.46 mag per 100 days. During the 2014 observing season, the 
decline rate is even slower, at 0.27 mag per 100 days in V and 0.30 
mag per 100 days in R. 

The MIR photometry available for SN 2009ip is quite limited, 
as can be seen in Fig. [5] From 2013 February until 2013 August, 
the decline in the 3.6 and 4.5pm bands appears consistent with that 
seen in the NIR, and only slightly slower than that seen at opti¬ 
cal wavelengths. Over the period 2013 July 25 to August 27, three 
epochs of MIR photometry were obtained; over this short period 
we observe a decline at 3.6pm, but an apparently constant flux at 
4.5 pm. It is possible that this is caused by contamination of the 
photometry of SN 2009ip by flux from a nearby source; as can be 
seen in Fig. [2] there are two bright sources close to SN 2009ip 
which are much brighter in the MIR than at optical wavelengths. 
Alternatively, if the constant 4.5 pm flux is due to SN 2009ip, then 
it may be caused by the formation of CO, which has its fundamental 
emission in this band. 

The colour of evolution of SN 2009ip is shown in Fig. [6] We 
find that in V-R, SN 2009ip has a roughly constant colour over the 
period from 300 to 800 days after the start of the 2012a event. In 
contrast, the B-V colour becomes slightly bluer over time. This is 
likely caused by the presence of emission from Fe group lines in the 
spectral region covered by the B-band filter. The V-H colour of SN 
2009ip can be measured until ~400 days, after this SN 2009ip is 
not detected in H-band, and we can only set upper limits to the NIR 
colour. The fact that SN 2009ip does not appear to become signifi¬ 
cantly redder after ~400 days would indicate that dust formation is 
not significant at this phase. 

We constructed a pseudo-bolometric BVRIJHKs lightcurve for 
SN 2009ip using the same technique as in iFraser et ail ( l2013ah . 
which is shown in Fig. [7] As expected given the steady decline of 
the individual photometric bands, the bolometric lightcurve shows 
a smooth decline. A comparison to the well-studied Type IIP SN 
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Figure 3. Sequence of optical spectra of SN 2009ip taken between 2013 April and December, with phases relative to MJD 56193 in parentheses. All spectra 
are shown in the rest frame, the wavelengths of key features are marked, while the positions of the strong telluric absorptions are shown with a © symbol. The 
WiFeS spectra and the EFOSC2 spectrum from November 25 are shown in light grey, with a binned version of the spectrum overplotted in black. A 4800K 
blackbody fit to the 2013 April 3 spectrum is shown in red. The higher resolution spectra taken during this period are shown in more detail in Fig. ED 
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8 M. Fraser et al. 


Table 2. Calibrated NIR magnitudes and MIR fluxes for SN 2009ip. 


Phase (d) 

Date 

MJD 

J (mag) 

H (mag) 

K (mag) 

3.6pm (mJy) 

4.5 p.m (mJy) 

Instrument 

130 

2013-01-31.01 

56323.01 

- 

- 

- 

0.1331 (0.0013) 

0.1628 (0.0012) 

IRAC 

201 

2013-04-11.37 

56394.37 

17.94 (0.06) 

17.92 (0.16) 

17.83 (0.11) 

- 

- 

SOFI 

303 

2013-07-23.16 

56496.16 

18.84 (0.18) 

18.99 (0.20) 

18.12(0.17) 

- 

- 

NOTCam 

306 

2013-07-25.95 

56498.95 

- 

- 

- 

0.0311 (0.0012) 

0.0460 (0.0009) 

IRAC 

326 

2013-08-14.55 

56518.55 

- 

- 

- 

0.0286 (0.0011) 

0.0464 (0.0008) 

IRAC 

339 

2013-08-27.74 

56531.74 

- 

- 

- 

0.0255 (0.0011) 

0.0474 (0.0008) 

IRAC 

341 

2013-08-30.02 

56534.02 

19.05 (0.16) 

19.41 (0.26) 

18.53 (0.24) 

- 

- 

NOTCam 

364 

2013-09-21.98 

56556.98 

19.44 (0.19) 

19.36 (0.26) 

<18.41 

- 

- 

NOTCam 

376 

2013-10-03.22 

56569.22 

19.18 (0.12) 

- 

- 

- 

- 

SOFI 

428 

2013-11-24.77 

56620.77 

19.68 (0.23) 

<18.84 

- 

- 

- 

NOTCam 

711 

2014-09-03.07 

56904.07 

20.67 (0.29) 

<20.5 

<18.8 

- 

- 

SOFI 

730 

2014-09-22.09 

56923.09 

20.78 (0.21) 

- 

- 

- 

- 

SOFI 



Rest wavelength (A) 

Figure 4. Sequence of optical spectra of SN 2009ip taken between 2014 April and September, with phases relative to MJD 56193 in parentheses. All spectra 
are shown in the rest frame, the wavelengths of key features are marked, while the position of the strong telluric absorptions are shown with a © symbol. The 
EFOSC2 spectrum from Aug 17/19 is shown in light grey, with a binned version of the spectrum overplotted in black. For comparison, the last spectrum taken 
of SN 2009ip in 2013 December (i.e. the final spectrum in Fig.[3j is also shown. The higher resolution spectra taken during this period are shown in more 
detail in Fig. ED 


2004et shows that the decline rate of SN 2009ip between —1250 
and +500 days from the start of the 2012a event is slower than typ¬ 
ical for a core-collapse SN at a comparable epoch. The degree to 
which radioactive Co contributes to the energy source of SN 2009ip 
during this late phase is, however, uncertain. As discussed later, the 
presence of strong and narrow Htr at this epoch suggests that CSM 
interaction is ongoing, and so one may expect any Co-powered tail 
to be supplemented by additional flux from the continuing con¬ 


version of kinetic energy. The strongest constraint on the Ni mass 
comes from the lightcurve in late 2012 December; from a compar¬ 
ison of the bolometric light curves of SNe 2009ip and 2004et at 
this epoch, we confirm the upper limit to the ejecte d Ni mass o f the 
former of Mn/ <0.02 M Q which was set bv lFraser et ali d2013al) . 

We also calculated bolometric light curves using the optical 
(BVRI), NIR (JHK), and MIR (LM) bands separately. The Spitzer 
3.6 and 4.5 pm fluxes were calibrated to the L and M bands us- 
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Table 3. Log of spectroscopic observations of SN 2009ip. 


Phase (d) 

Date 

MJD^ 

Telescope 

Instrument 

Grism 

Range (A) 

193 

2013-04-04.39 

56386.39 

NTT 

EFOSC2 

Grl3 

3650-9246 

201 

2013-04-12.00 

56394.00 

ANU 2.3m 

WiFeS 

B3000+R7000 

4171-7021 

202 

2013-04-13.40 

56395.40 

NTT 

EFOSC2 

Grll 

3346-7466 

203 

2013-04-14.39 

56396.39 

NTT 

EFOSC2 

Grl 6 

5995-9990 

207 

2013-04-18.39 

56400.39 

NTT 

SOFI 

BG 

9277-16345 

209 

2013-04-20.39 

56402.39 

NTT 

EFOSC2 

Grll 

3347-7467 

222 

2013-05-03.40 

56415.40 

Gemini-S 

CMOS 

R400 

4208-8439 

226 

2013-05-06.35 

56418.35 

VLT 

XShooter 

UVB+VIS+NIR 

3000-24800 

239 

2013-05-20.35 

56432.35 

Gemini-S 

CMOS 

R400 

4207-8439 

242 

2013-05-23.43 

56435.43 

Gemini-S 

CMOS 

B600 

3191-6076 

245 

2013-05-26.00 

56438.00 

ANU 2.3m 

WiFeS 

B3000+R7000 

3500-7023 

264 

2013-06-14.36 

56457.36 

Gemini-S 

CMOS 

R400 

4208-8437 

272 

2013-06-22.00 

56465.00 

ANU 2.3m 

WiFeS 

B3000+R7000 

4169-7023 

280 

2013-06-30.43 

56473.43 

Gemini-S 

CMOS 

R400 

4210-8441 

303 

2013-07-23.38 

56496.38 

Gemini-S 

CMOS 

R400 

4207-8440 

323 

2013-08-12.19 

56516.19 

Gemini-S 

CMOS 

R400 

4208-8439 

353 

2013-09-12.05 

56547.05 

NTT 

EFOSC2 

Grl 1+Grl6 

3343-9986 

374 

2013-10-03.14 

56568.14 

VLT 

FORS2 

300V 

4226-9639 

428 

2013-11-26.08 

56622.08 

NTT 

EFOSC2 

Grll 

3365-7478 

434 

2013-12-02.06 

56628.06 

VLT 

FORS2 

300V 

4224-9632 

575 

2014-04-22.37 

56769.37 

NTT 

EFOSC2 

Grl 3 

3650-9246 

589 

2014-05-06.36 

56783.36 

NTT 

EFOSC2 

Grll 

3365-7478 

590 

2014-05-07.33 

56784.33 

NTT 

EFOSC2 

Grl6 

5995-9990 

590 

2014-05-07.39 

56784.39 

VLT 

XShooter 

UVB+VIS+NIR 

3000-24800 

693/5 

2014-08-17/19 

56887.16 

VLT 

XShooter 

UVB+VIS+NIR 

3000-24800 

726 

2014-09-20.36 

56920.36 

Keck 

DEIMOS 

600ZD 

4200-7640 


ing the zeropoints in lCampins et al.l l ll985h . For the period where 
there were MIR observations (2013 January - August), the contri¬ 
bution at these wavelengths to the total optical-MIR luminosity of 
SN 2009ip was only ~10 per cent. The NIR contributed a further 
~30-35 per cent of the luminosity, while the remainder of the flux 
emerged at optical wavelengths as shown in Fig. [8] We see that af¬ 
ter 2013 September, the fraction of flux in the NIR relative to the 
optical appears to decrease by about 5 per cent. It is possible that 
this is due to a large amount of the flux in SN 2009ip emerging in 
line emission, much of which is at optical wavelengths. 

iMartin et~aD d2013l) and lMargutti et al.l J20 1 4) put forward the 
intriguing suggestion that there was a periodic component of the 
lightcurve of SN 2009ip. With a view to testing this in the late 
time data, we took the R-band lightcurve of SN 2009ip between 
2013 April 3 and December 24, and made a least-squares linear fit 
to the data. To remove the overall decline from the lightcurve, a 
“residual lightcurve” was constructed, consisting of the measured 
magnitude at each epoch minus the interpolated magnitude from 
the linear fit. The R-band data were used for this as this was the 
filter which had the most data points (36 epochs). We then used the 
Lomb-Scargle algorithm dLombl 1 976: Scar glel 19821) to derive a pe- 
riodogram from the residual lightcurve. The resulting periodogram 
is shown in Fig. [9] but as expected from visual inspection of the 
lightcurve, we see no significant peaks in the power spectrum that 
would suggest a periodic signal. The absence of any periodic com¬ 
ponent to the R-band light curve may indicate that the region of 
CSM which the ejecta is now interacting with has a different den¬ 
sity profile, or was formed through a different mass loss process 
to that which formed the CSM closer to SN 2009ip. The progen¬ 
itor candidate identified by iFolev et all d201 lh for SN 2009ip had 
a magnitude of 21.84±0.17 mag in F606lV-band images taken in 
1999 June. On 2014 December 23, the measured magnitudes of SN 


2009ip were V = 21.73 ± 0.17 andR=20.55 ± 0.07 mag. In V-band, 
SN 2009ip reached a similar magnitude to the pur ported proge nitor 
in late 2014. However, as the HST images used bv lFolev et alj were 
taken with the F606W filter, they include the wavelength of Hff. As 
the H a emission seen in SN 2009ip is strong, even during quies¬ 
cent periods l lPastorello et al.ll2013l) . it is likely that Hff emission is 
contributing to the flux here also. Hence the most appropriate filter 
to make a comparison to the pre-discovery HST images is R-band, 
which similarly includes Hff. SN 2009ip is still ~1.3 mag brighter 
than the progenitor candidate in R-band, and assuming a constant 
decline rate of ~0.3 mag per 100 days, will only reach a comparable 
magnitude around 2016 March. 

4 SPECTROSCOPIC EVOLUTION 
4.1 Optical spectra 

As is evident from Figs. [3 and [4] SN 2009ip does not undergo any 
dramatic spectral evolution between 2013 April and 2014 Septem¬ 
ber. From the first spectrum of 2013, which was taken on April 3, 
until June 22, we do not observe the emergence of any new lines 
or significant changes in velocity or line profile shape. In the fol¬ 
lowing, we describe the 2013 April 3 spectrum as representative of 
SN 2009ip over the subsequent three months. We note that the five 
highest resolution and signal-to-noise spectra are shown in more 
detail in Fig. m All of the reduced spectra presented in this pa- 
per w ill be placed in the the WiseRep archival dYaron & Gal-Yaml 

s. 


5 http ://www. weizmann. ac. il/astrophysic s/wiserep/ 
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Sep 2012 Dec 2012 Mar 2013 Jun2013 Sep 2013 Dec 2013 Mar 2014 Jun2014 Sep 2014 Dec 2014 



Figure 5. Optical and infrared lightcurves for SN 2009ip. The data from lPastorello et ~all ^20 1 3|) and lFraser et aljf2013al) covering the 2012a and 2012b events 
are shown with open points. The black line shows the typical 0.98 mag per 100 days decline expected during the radioactively powered tail phase of a core¬ 
collapse SN. Upper limits where SN 2009ip was not detected are indicated with arrows; the green horizontal arrow at the lower right indicates the presumed 
quiescent F606W (~R-band) magnitude measured from archival HST images by I Foley et ali 1201 1.1. scaled by -2 mag to match the R-band lightcurve. The 
source of each point in the optical lightcurve is indicated by the key. 

As shown in Fig. [3] the spectrum from April 3 has a rela¬ 
tively flat continuum. The most prominent feature in the spectrum 
is H<r emission, which is well fitted with a Lorentzian profile of 
FWHM 1450±100 kms -1 . Fla- appears to be slightly asymmetric, 
as shown in Fig. QT| with a deficit of flux in the blue wing. We 
also see H/3 with a similar velocity and line profile, and a weak Hy 
line. While the line profile of Flo- appears to narrow slowly over 
the course of our monitoring campaign, for the first three months 
this evolution is hard to discern. The Balmer decrement (L Ha /L H p) 
can be used to constrain the physical environment of SN 2009ip 
dLevesque et alJl2014l ; lFraser et al.ll2013al) . and in the April spec¬ 
trum of SN 2009ip we measure a Balmer decrement of ~12. This 
is much higher than in 2012, and could indicate Case C recombi¬ 
nation, where Ha- is optically thick (IXu et al.ll 1992h . However, the 
lack of any clear absorption in the Balmer lines would disfavour 
this scenario. Alternatively, if Har is optically thin, a lower elec¬ 
tron temperature and a higher density than at peak^could perhaps 
give rise to the observed decrement forake & UlrichlH 98(j F1. The 


6 In lFraser et al.l i2013al) . we did not use the results of lDrake & UlrictJ . as 
the high luminosity of the 2012b event requires the inclusion of the radia¬ 
tion field when computing line fluxes. Here, as SN 2009ip is much fainter, 


higher density could be explained as the photosphere recedes in¬ 
wards through material ejected in a wind or mass loss with a oc R 2 
dependence, although it is unclear whether such a simple scenario 
would necessarily be appropriate to SN 2009ip. The strong Her line 
is also interesting in the context of the low photospheric temper¬ 
ature. As the spectral energy distribution of SN 2009ip shifts to 
cooler temperatures than at peak, there are no longer the UV pho¬ 
tons necessary to produce recombination lines. It is likely that X- 
rays from shocks in the interaction region are now supplying the 
necessary energy to excite Her, and perhaps also help maintain the 
relatively constant blue colour of SN 2009ip (Fig. In the red 
part of the spectrum, we identify some of the higher order Paschen 
lines. These can be seen more clearly in the spectrum from April 
13, where Pa lines as high as the 12—>3 transition at /18750 can be 
seen close to the Ca 11 NIR triplet. 

To estimate the temperature of SN 2009ip in 2013 April, we 
fitted a blackbody function to the spectra obtained in 2013 April, 
with the exception of the WiFeS data. The regions of the spec¬ 
trum around Her and H/J were excluded from the fit in all cases. We 


the assumption of lDrake & Ulrich! that the radiation field is weak is more 
appropriate. 
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Figure 10. Rest frame high resolution VLT+XShooter, VLT+FORS and Gemini+GMOS spectra of SN 2009ip, taken between 2013 May and 2014 May. In the 
top panel, the VLT+XShooter spectrum from 2014 May 7 has been multiplied by 3 for clarity, all other spectra have been offset in flux as necessary. The 2013 
May 6 and 2014 May 7 XShooter spectra have also been boxcar smoothed using a 15 pixel window (red and black lines) to emphasise the weaker features. 
BO^Q^^nd the region of atmospheric absorption between 8900 and 9900 A are marked with a © symbol. 
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Figure 6. The B-V, V-R and V—H colour evolution of SN 2009ip between 
2013 April and 2014 December. Limits on the colour of SN 2009ip, from 
epochs where it was only detected in a single filter, are indicated with ar¬ 
rows. The colour curves have been corrected for foreground extinction. 
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Figure 7. BVRIJHK p seudo-bolometric lightcu rve of SN 2009ip (open 
points; IPastorello et al Jl 201 3t iFraser et alj|2013al) , plus the additional late¬ 
time data presented in this paper ( solid points). For com parison we 
show the core-collap se SNe 2005cs jPastorello et alj [20091) and 2004et 
( Maguire et al. 2010) shifted so that their explosion epoch matches the start 
of the 2012a event. 
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Figure 8. The fraction of flux in optical ( BVRI ), NIR ( JHK) and MIR 
{Spitzer) bands for SN 2009ip during the tail phase. Epochs later than 2013 
December have scant NIR coverage (SN 2009ip is detected in 7-band only) 
and hence are not shown; the time period covered by Spitzer MIR observa¬ 
tions is also indicated. 
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Figure 9. Lower panel: R-band lightcurve of SN 2009ip between 2013 April 
and December, after the subtraction of a linear least-squares fit to the data 
points. Upper panel: Lomb-Scargle periodogram of the same, where cj is 
the angular frequency in radians per day, and power is in arbitrary units. 
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Figure 11. Evolution of the Ho line profile over time. The sequence con¬ 
tains all optical spectra from Table [3] which cover the wavelength of Ho, 
with the exception of the EFOSC2 Gr#13 spectra which were of lower res¬ 
olution. The date of each spectrum is indicated, as is the phase relative to 
the start of the 2012b event, MJD=56193. All spectra have been convolved 
with a Gaussian to degrade them to a common resolution of R=1000, have 
been normalised to the peak of Ho-, and have been shifted so that the peak 
of Her is at 0 A. The emergence of a red shoulder in Ha at at the end of 2014 
is apparent in the final two spectra. 


Table 4. Measurements of the FWHM of Ha in SN 2009ip, as found from 
a Gaussian fit to the emission. 


Phase (d) 

Date 

FWHM (A) 

202 

2013-04-12 

34.3 

201 

2013-04-13 

37.1 

203 

2013-04-14 

36.3 

209 

2013-04-20 

36.7 

222 

2013-05-03 

35.1 

226 

2013-05-06 

33.1 

239 

2013-05-20 

34.4 

245 

2013-05-26 

32.7 

264 

2013-06-14 

33.2 

272 

2013-06-22 

32.3 

280 

2013-06-30 

33.0 

303 

2013-07-23 

33.0 

323 

2013-08-12 

32.7 

374 

2013-10-03 

32.2 

428 

2013-11-26 

31.2 

434 

2013-12-02 

30.1 

589 

2014-05-06 

28.2 

590 

2014-05-07 

28.7 

590 

2014-05-07 

26.2 

726 

2014-09-20 

25.3 


measure an effective temperature of 4800 K at this epoch (—h 200 
d). Subsequent spectra show a slight increase in flux to the blue 
of ~5500 A, but this is likely caused by a relative increase in the 
strength of the Fe n multiplets between 5000 and 5500 A. We ob¬ 
serve a similar effect in the B-V colour. At later epochs, a clear 
continuum arising from a photosphere cannot be discerned, and so 
we cannot measure a physically meaningful effective temperature. 


In addition to H, we see He, Na, Ca, and Fe in the spectrum of 
SN 2009ip in 2013 April. Weak He i 47065 is also present in emis¬ 
sion. The other transitions of He in the optical are unfortunately 
blended with other lines; He i 45876 is blended with Na i D, while 
both the 44922 and 45015 He i lines lie in the region of strong 
emission from Multiplet 42 of Fe. We alsjj_tentatively detect He i 
46678 which was seen at earlier epochs dFraser et al.ll2013ah . The 
He i 47065 line is not particularly clear in the April 3 spectrum, 
and so we instead measure its FWHM in the April 14 spectrum 
taken with EFOSC2 + Gr#16, where we find it to have a FWHM 
of only ~1000 kmr 1 , which is significantly slower than Ha. The 
0 1 48446 line can also be seen to emerge in the April 13 spectrum, 
with a FWHM of ~1250 km s _1 . This feature was possibly present 
in the earliest spectrum taken on April 3, albeit blended with the 
blue wing of Ca ii 48498. There is perhaps some contribution to the 
line from the 18—>3 transition of H (48438), although the intensity 
of the feature in comparison to the other Pa lines suggests that H is 
not the dominant element contributing to any blend. Furthermore, 
the feature grows in strength over the following six months, consis¬ 
tent with the line being formed primarily by O. O 1 48446 emission 
in Type II SNe typically arises via Lyman j3 pumping in O i 41026 
(e.g. m Lvall 19931 ). The strength of this line in Type II SNe has been 
interpreted as radiative cros s-talk between H and O zones dOlival 
119931 : Ijerkstrand et al.ll2012h . although unless it is very strong, the 
48446 line (and O i 1.129pm) can be produced Ly-beta scattering 
in primordial material. 

The only line with a clear absorption component is Na i D, 
which has a P Cygni profile with a minimum at -3000 km s _1 . The 
peak of the emission is at 5896 A, indicating that the emission is 
definitely caused by Na, although we caution that the absorption 
may be contaminated by He i 45876. There is also a tenuous de¬ 
tection of a second absorption component in Na at ~15000 km s -1 , 
which would match to the two component absorption seen in ear¬ 
lier epochs dFraser et al.ll2013ah . Na i D tends to remain optically 
thick for a long time, as it is a resonance line from the ground 
state, and so it is unsurprising that we continue to see absorp¬ 
tion. [Ca ii] 447291,7324 is present, as is the Ca n NIR triplet in 
pure emission. [Ca ii] appears broader than other features, which 
is due to the blending of the two components in a doublet, and 
possibly due to some contamination from He i 47281. We mea¬ 
sure a FWHM of 1700 km s _1 for the 48662 line in the Ca ii NIR 
triplet, similar to the H and He velocities. Aside from these lines, 
the only other features present are from Fe, where we see several 
permitte d multiplets in Fe n, the strongest of which is Multiplet 42 
dMoordll945l) ; although there is undoubtedly contamination from 
He i 44922 and 45015. We also see a possible weak feature at 
~7720 A, although this region is affected by telluric absorption. 
The feature is offset by -2000 km s _1 from the rest wavelength of 
O i 47774, and an alternative identification could be Fe i 47710. 
The permitted Fe and Na, and the permitted and forbidden Ca can 
all likel y be accounted for by primordial metals in the envelope of 
the star dLi & McCravll 1993 k iMaguire et al.ll2012l ; Ijerkstrand et al.l 
EMU we s tress that these elemen ts were all seen at comparable 
velocities by lPastorello et all d2013l) in spectra of SN 2009ip taken 
two years prior to the 2012a and 2012b events, and hence cannot 
be unambiguously associated with the products of SN nucleosyn- 


7 Although Magui re"et al J d2012l) and Ijerkstrand et alj 120121) presented 
models for the lower mass red supergiant progenitor of a Type IIP SN, the 
conclusion that much of the line emission seen at late times arises from the 
envelope should also apply here. 
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thesis. The similar velocities measured for H, He and Ca are also 
consistent with these lines forming in the same location. We also 
note that whatever the identification of the ~7720 A feature, both 
Fe 1/17710 and O i 47774 were seen in the September 2009 spectra 
JPastorello et al.lEoi3l) . The only ambiguous feature in the spectra 
of SN 2009ip in 2013 April is a line around ~4580 A, which could 
possibly be identified with Mg i] 44571. The feature is centered 
600 kms~’ to the red of the rest wavelength of Mg. As discussed 
later, we do not see the Mg i 1.5 pm line in the NIR spectra taken 
around this epoch, although the NIR line is a recombination rather 
than a cooling line. 

In the XShooter and GMOS spectra obtained in 2013 May, 
we observe the first evolution in the spectra of SN 2009ip, with 
the strengthening of the 48446 O i line (Fig. 110b. which was only 
marginally detected in 2013 April. By 2013 September, the line 
has grown in strength to be comparable with the components of the 
close by Ca n NIR triplet. In the June 30 spectrum, we see the emer¬ 
gence of [Fe n] T7155. The feature is quite weak, but continues to 
grow in strength until our latest spectrum in December. Interest¬ 
ingly, the 47155 feature was not seen i n the spe ctra covering the 
pre-2012 evolution of SN 2009ip dPastorello et al.ll2013l ). while the 
line is often seen during the nebular phase of core-collapse SNe. 

The final significant change in the spectrum of SN 2009ip 
in late 2013 is the emergence of what appears to be weak [O i] 
446300,6364 in the VLT+FORS spectrum taken on October 2. 
Emission lines consistent with both components of the doublet are 
seen, an appear to have approxim ately the same intensity. These 
lines were also seen by IFox et all 120151) in their spectrum of SN 
2009ip from 2014 June. The lines are quite weak but clearly real; 
they are detected in two separate spectra, most clearly on Octo¬ 
ber 2, where the spectrum has a S/N in the continuum of ~25. The 
lines are also comparable in strength to the nearby Multiplet 74 
of Fe n, as can be seen in Fig. ED As an alternative identification, 
the presumed O lines could be associated with Fe n, however we 
were unable to find a consistent match to any of the common mul- 
tiplets dMoorell 1 9451) . Multiplet 40 has a line at 6369 A which is a 
reasonable match, however there should also be stronger emission 
from the same multiplet at 46432 (which is possibly blended with 
the line at 6445 A), and at 46516 (which may be lost in the blue 
wing of Hu. Multiplet 200 has a component at 46305, but again 
this should be accompanied by a line of approximately the same 
strength at 46175 which we do not observe. In the October 2 spec¬ 
trum, we fitted Lorentzian profiles to the two components of the 
O doublet and measured an average FWHM of 1400 km s _1 while 
in the December 1 spectrum we measured 1200 kms -1 and 1400 
km s _1 for the 46300 and 46364 lines respectively. This velocity is 
similar to that seen in the nearby Fe Multiplet 74 at 1400 kms -1 , 
and consistent with that of Her. 

A ratio of 1:1 between the intensity of the two components of 
[O i] 446300,6364 would suggest optically thick conditions where 
the lines are forming; if the lines are forming in the ejecta, then 
the line ratio would be expected to increase to 3:1 over time as the 
ejecta expands and the density decreases (e.g lchugailll992l) . Pri¬ 
mordial [O i] 446300,6364 is always observed in the 3:1 ratio, and 
hence an observed ratio of ~ 1:1 is suggestive of synthesized oxy¬ 
gen which has a higher density. However, as can be seen in Fig. 1 121 
its also quite possible that there is contamination of the 46364 line 
by the many iron-group lines in this region. Furthermore, in core¬ 
collapse SN spectra the [O i] lines are much stronger than these 
surrounding Fe lines; and so if this is synthesized oxygen its mass 
must be very small. In later subsequent spectra taken in 2014, these 
lines are not detected, as can be seen in Fig. ED 
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Figure 12. The region of the presumed [O i] 446300,6364 doublet in two 
late time spectra of SN 2009ip, taken in 2013 October and December. The 
rest wavelengths of the two components of the doublet, along with those 
of the nearby Fe n Multiplet 74 are indicated. The 2014 May 7 spectrum 
is also shown in light blue, but even after binning (dark blue line) the [O i] 
doublet is not visible at this epoch. 


The XShooter and GMOS spectra obtained in 2013 October 
also reveal some structure in the Ca lines. As can be seen in Fig.1101 
all three lines of the Ca II NIR triplet show an additional emission 
component, offset by ~550 km s _1 to the red. This feature is real, as 
it can be seen in multiple spectra at different epochs and from dif¬ 
ferent instruments. It is unclear whether the redshifted components 
are also present in the [Ca n] 447291,7324 lines; as these lines are 
only separated by ~30 A, they appear blended together and any 
additional component would be hard to disentangle. A similar red- 
shifted emission component is not seen in any of the other elements, 
with H, He and Fe all showing relatively symmetric profiles. The 
lines of O are weak, while Na i D is blended with He, making it dif¬ 
ficult to assess their profiles, however neither of these lines shows 
an obvious emission shoulder either. 

SN 2009ip was behind the Sun during the period from 2014 
January to March. When SN 2009ip was first observed again in 
2014 April, it appeared to be relatively unchanged from the previ¬ 
ous spectrum in 2013 December (Fig. 0. Interestingly, and con¬ 
trary to naive expectations if SN 2009ip had undergone core¬ 
collapse, the [O i] 446300,6364 feature had not strengthened since 
its first appearance in 2013 October. In fact, only the 46300 com¬ 
ponent of the doublet was detected in the XShooter spectrum from 
2014 May. Whether the naive expectation that [O i] would continue 
to strengthen is, of course, uncertain; the models of ljerkstrand et al.l 
J2014t) only cover the period up to ~450 days after core-collapse, 
by which point the [O i] lines appear to have reached their maxi¬ 
mum strength. 

The subsequent spectrum we obtained of SN 2009ip using 
VLT+XShooter in 2014 May also revealed an additional emission 
component in the red wing of the Balnrer lines. The feature is seen 
to be offset by +500 km s' 1 from the rest wavelengt h of H o-, and is 
also seen as a faint shoulder in H/J. IGraham et alj 1 120141) also saw 
a slight excess of flux in the red wing of Hu in as early as 2013 
June, although at this stage the asymmetry was less clear than the 
distinct shoulder seen in our data from 2014 May. The velocity of 
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the red emission feature in Ha is similar to that seen in the sec¬ 
ond component of the Ca II NIR triplet in 2014 October; we note 
that the redshifted Ca emission is also still present at this epoch. 
We also note that the He emission lines in the 2014 May XShooter 
spectrum are extremely narrow, with He i 47065 having a velocity 
of only ~ 180 km s _1 as measured from the FWHM of a Lorentzian 
fit to the line. This is marked contrast to the velocity of 750 km s _1 
measured in the same line in the XShooter spectrum taken a year 
earlier in 2014 May. It is also significantly lower than the velocities 
seen in the Balmer lines, which have FWHM of around 1000 km s -1 
in 2014. We also see similarly low velocities in the line at 5016 A 
suggesting that this is dominated by He 1 45016 at this stage, rather 
than by Fe n Multiplet 42. 

Our final spectrum of SN 2009ip was obtained on 2014 
September 19, and is similar to those from 2014 April and May. 
The spectrum has an excess of flux in the blue, which we attribute 
to a pseudo-continuum formed of many weak Fe-group lines. Ha 
and H ft are still present, with narrow (~1000 kms~') Lorentzian 
profiles. [Ca ii] is present, while it is unclear whether Na is still 
detected due to the He emission at 5876 A which could mask the 
presence of Na i D. Finally, the He lines are extremely narrow, with 
velocities of only ~200 km s -1 . The lower velocities for the He lines 
when compared to H suggest that these may arise from different re¬ 
gions. Such relative strength of the He lines may also point to a 
CSM which is richer in He than a typical Type Iln SN. 

We measured the FWHM and central wavelength of Ha in all 
the spectra obtained of SN 2009ip; these measurements are listed in 
Table[4] Each spectrum was convolved to a common resolution, and 
a Gaussian was fit to Ha, allowing the central wavelength, FWHM 
and a normalisation constant to vary as free parameters. We see no 
systematic shift in the position of Ha over time, suggesting that 
substantial dust formation is not occurring. The measured values 
for the FWHM of Ha are plotted in Fig. |T3] From Fig. |T3] the 
decrease in FWHM of Ha over time, from ~1600 kms~' in 2013 
April, to — 1100 km s_[ in 2014 September, is clear. 

While ISmith et ahl d2013l) suggested that dust formed in the 
material lost by SN 2009ip during previous eruptions, we find no 
evidence for dust formation since the peak of the 2012b event. In 
the spectra obtained for SN 2009ip between 2013 April and 2014 
September, Ha appears to be centred on its rest wavelength, with 
no systematic blueshift which would be indicative of dust. Further¬ 
more, we see no sign of a deficit of flux in the red wing of Ha. 
We also do not observe any shifts or deficit in the red wing of any 
of the other lines. The H-K colour is positive, while J-H is nega¬ 
tive, and the V-K colour is only 1.5, providing ambiguous evidence 
for a K-band excess. Furthermore, as the flux in the continuum of 
SN 2009ip is relatively weak compared to the strong line emission, 
it is possible that broad band photometry does not provide a good 
estimate of the shape of the underlying SED. 

4.2 NIR spectra 

The NIR spectra of SN 2009ip are shown in Fig. |T4] while in Fig. 
I15lwe show a zoomed in region of the 2013 May XShooter spec¬ 
trum, covering the location of the first overtone of CO. The spec¬ 
trum taken with the NIR arm of XShooter in 2014 May suffers from 
extremely low S/N, and hence is not discussed here. The two spec¬ 
tra from 2013 were taken less than a month apart (with no apparent 
evolution between them), and are dominated by H, and to a lesser 
extent He. We see the Pa lines, including Paa which is of sufficient 
strength to be recovered even with the strong high atmospheric ab¬ 
sorption between 1.8 and 2.0 pm. Lines from the Brackett series are 
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Figure 13. The central wavelength and FWHM of Ha. These values were 
measured for all spectra via a Gaussian fit to the line. 

visible in the XShooter spectrum. We also observe the two strongest 
He i lines at 1.08 and 2.05 pm. 

The Mg i line at 1.5 pm is not observed in either of the NIR 
spectra, although we note that the S/N per pixel in the continuum 
is ~3 for the 2013 May 6 spectrum. We take the measured flux of 
2.5xl(L 17 erg s~' errT 2 in the weakly detected (and nearby) 14—>4 
Brackett line as a conservative upper limit to the flux in the Mg 
line. The measured flux of the possible Mg i] 44571 line at this 
epoch is 1.6xlCT 15 erg s _1 cnT 2 implying that the latter is at least 
~60 times stronger. This implies that that the formation of Mg i] 
44571 and Mg 1 1.5 pm are different, with the former a cooling line 
from primordial gas, while the latter is a recombination line with 
no strong contribution from primordial Mg. 

The extreme red part of the XShooter spectrum shows an in¬ 
crease in flux at around 2.25pm and beyond. Fig.[l5]shows a com¬ 
parison of this to the Pfund series in the limit of 5—» oo, also shown 
are the bandheads from molecular CO emission. The observed in¬ 
crease in flux is roughly consistent with the position of these, and 
while the flux calibration of the observational data is least reliable 
at the edges of spectra, it seems plausible that CO has begun to 
form by ~280 days after the start of the 2012a event. If this is the 
case, then the measured 4.5 pm fluxes from Spitzer + IRAC pre¬ 
sented in Section [2] should be regarded as upper limits to the con¬ 
tinuum flux, as the filter covers the wavelength of the fundamental 
emission from CO. We note that_CO band heads have been seen 
in before in cool, evolved stars l lBieging et alj2002h : while CO has 
also been seen in Eta Car, where it is b elieved to have formed in 
CNO-processed material iLoinard et aljhoil) . From this, the pos¬ 
sible presence of CO emission here cannot be taken as evidence 
either for or against a core-collapse. 


5 DISCUSSION AND CONCLUSIONS 

Despite the fact that SN 2009ip is one of the best monitored ex- 
tragalactic transients to date, even five years after its initial discov¬ 
ery it continues to present a puzzle. Perhaps the most significant 
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Figure 14. NIR spectra of SN 2009ip. To aid in the identification of lines, a binned version of each spectrum is shown in black, while the unbinned spectra are 
in grey. The red shaded areas correspond to low atmospheric transmission. The final XShooter spectrum taken in 2014 May contains essentially no flux in the 
NIR from SN 2009ip, with the possible exception of some line emission from He 1.08 pm. 



Rest wavelength (pm) 

Figure 15. A'-band region of the XShooter spectrum from 2013 May 6 (in 
grey) showing the increasing flux at red wavelengths, along with a smoth¬ 
ered, binned version of the same spectrum shown with a black line. The 
positions of the CO 3 —* 2 and 2 —* 0 band heads are also shown, along 
with the Pfund series and Bry and He i 2.06pm. 


among these questions is whether SN 2009ip has indeed exploded 
as a core-collapse SN. 


5.1 The energetics of SN 2009ip 

The lightcurve of the 2012a event appears to follow a steady rise 
to peak, followed by a symmetric decline. This is reminiscent of 
the diffus ion peak seen in Type I SNe. In fact, in the scenario pro¬ 
posed by iMauerhan et ]d] d2013h for SN 2009ip, the 2012a event 
was the weak explosion of SN 2009ip as a core-collapse SN (while 
the 2012b event was caused by the ejecta from this event hitting the 
CSM). 

Following l Arnetli d 199fil) . the diffusion timescale for a homolo¬ 
gously expanding, spherical body of uniform density can be shown 
to be approximately 


U - 


/ 3ftM ej \ 1/2 
\47TCVej/ 


( 1 ) 


where M s] and v e j are the ejecta mass and velocity respectively, k is 
the opacity of the ejecta, and c is the speed of light. If we take the 
~20 day rise of the 2012a event as being characteristic of the diffu¬ 
sion time, and a velocity of 8500 km s~' (from the mini ma o f the P 
Cygni absorptions in the Balmer lines: IPastorello et al.ll2013i ). and 
assume an opacity /c=0.35 cm 2 g -1 (assuming e scat tering dom¬ 
inates in the H-rich ejecta; iRvbicki & Lightmanll 19791 ), the ejecta 
mass of the 2012a event is only ~0.44 M 0 . 

We note that there are two caveats to this estimate, namely 
whether such a simplistic estimate is appropriate here, and whether 
the absence of photometric coverage in the UV bands during the 
rise of 2012a event has a significant influence on the bolometric 
lightcurve. The optical flux after maximum in the 2012a event was 
a factor ~6 greater than the UV flux JFraser et al.l20133) . and while 
the effective temperature prior to the 2012a event maximum was 
higher, it did not reach the 17000 K seen in the 2012b event when 
the UV flux contribution did dominate prior to maximum light 
dFraser et alj) . All in all, we consider it unlikely that the shape of 
the 2012a event lightcurve would be drastically different with wider 
wavelength coverage. Whether such a simplistic estimate of the 
lightcurve is suitable in the case of SN 2009ip is less certain. In par¬ 
ticular, the Arnett model assumes a centrally located power source 
for a transient, while in the case of SN 2009ip, this may be supple¬ 
mented by additional energy from interaction with a CSM. Non- 
spherical geometry may also significantly affect the behaviour of 
SN 2009ip. However, and with these caveats in mind, a toy model 
may still provide a useful guide as to the energetics and masses in¬ 
volved. In particular, if the 2012a event is the core-collapse of SN 
2009ip, then the ejecta mass estimate should be reasonably valid, 
and at least should provide an upper limit. 

If we assume that 0.44 M 0 was ejected at ~8500 km s~* at the 
start of the 2012a event, and that the 2012b event was powered by 
the kinetic energy of this material, then the ~ 50 days until the start 
of the 2012b event would imply the inner edge of the CSM lies at 
4x 10 15 cm (270 AU) from the progenitor. 

Toy models have been used to estimate the ejecta and 
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Date 


Sep 2012 Mar 2013 Sep 2013 Feb 2014 Aug 2014 



Figure 16. Bolometric lightcurve of SN 2009ip, compared to the energy 
from the decay of 0.015 M e of 56 Co (based on the upper limit to the mass 
of 56 Co ejected by SN 2009ip from lFraser et all 12013 2); assuming full trap¬ 
ping of y-rays in the ejecta, and to a power law fit to the late time lightcurve 
with exponent —1.74). 


CSM characteristics of int e racting SNe (e.g. Chu gai & Da nzigerl 
1994b Moriva et all 1201 3l : IChevalier & Franssonl I 1 994^ . From 

J <2 ' 


Moriva et alj d20 1 3l) . we see that in the case of steady mass loss 
(p oc R 2 ), the late time luminosity decline should follow oc r 1 ' 5 
(eq. 29). In the case of SN 2009ip, we measure a decline with a 
time dependence oc f -1 - 74 , as shown in Fig. |T6] The disagreement 
with the expected slope for steady mass loss is unsurprising, given 
the observed eruptive behaviour of SN 2009ip prior to the 2012a 
and 2012b events. 

The emergence of the red shoulder seen in the Her line at late 
times (shown in more detail in Fig. 1171 is evidence that there is still 
structure in the CSM around SN 2009ip at this epoch. A possible 
formation mechanism for the feature could be interaction between 
ejecta and CSM on the far side of SN 2009ip, where the energy 
from X-rays created at the shock interface is absorbed locally; i.e. 
the emergence of a new kinematical component. In this case, the 
absence of a corresponding feature in the blue wing of Her indicates 
asymmetry about SN 2009ip along the line of sight (consistent with 
the findings of lMauerhan et al.l2014l from spectropolarimetry). We 
note that a similar line profile was seen for some Type Iln SNe, in¬ 
cluding SN 1995N dPastorello et alJl201ll) . The decrease in promi¬ 
nence of the wings of Her over time is indicative of reduced electron 
scattering. This is a consistent with the decreasing energy input, 
giving rise to lower ionization. 


5.2 Possible nucleosynthesis in SN 2009ip 

The presence of strong, broad O in late-time, nebular spectra 
is a classical signature of a core-collapse SN (CCSN). [O i] 
446300,6364 is formed from both primordial O in the envelope, 
and synthesised O from the core of the progenitor, and the doublet 
is seen to be stronger for more massive progenitors (which contain 
more massive cores, e.g. ljerkstrand et alj2014l) . The appearance of 
weak [O i] in the spectrum of SN 2009ip approximately a year af¬ 
ter the 2012a event could hence be seen as evidence that SN 2009ip 
did indeed become a genuine SN. However, if this were the case, 
then we would expect [O i] to become progressively more promi¬ 
nent over time, as the densit y of the line emitting region in the SN 
core decreases (e.g. ljerkstrand et al]|2012h . In contrast, we see that 



-50 -25 0 25 50 

Wavelength (A) 


Figure 17. The profile of Hi? as seen in the XShooter spectrum obtained on 
2014 May 7 (black), compared to the line profile observed for Hu one year 
earlier on 2013 May 6 (red). The wavelength scale is shown relative to the 
rest wavelength of the line, the red emission shoulder discussed in the text 
is clearly visible at +550 km s~ 1 . 


the [O i], and indeed the other lines which are plausibly associated 
with a CCSN such as [Fe n] 47155, remain weak. 

We also note that even without nucleosynthesied material 
from either a SN explosion or stellar evolution prior to collapse, 
radiative transfer models of nebular Type II SN spectra can show 
emission from C, Na and O solely arising from primordial material 
which was present in the stellar envelope. In Fig. |T8] we plot the 
emission from the inner H zone at 400 days, from the CCSN model 
of llerkstrand et ah ( 2012 ). This model shows a reasonable match 
to the observed 2013 October 2 spectrum of SN 2009ip, which was 
taken at a similar epoch relative to the start of the 2012a event. The 
[O i] /T16300,6364 lines are present in the ljerkstrand et aD d2012h 
model, and are actually stronger than those observed in SN 2009ip. 
In this model, cooling of the H gas occurs through metal lines such 
as [Ca n] 447291,7323, the Ca ii NIR triplet, [Fe n] 47155, [O i] 
446300,6364, and Mg i] 44571, while recombination gives rise to 
lines from H, He, and some of the metals. The presence of these 
lines emphasises the need for detailed modelling to distinguish a 
spectrum arising from H gas with metals at solar abundances (as 
in a wind or non-terminal outburst) and that of such gas plus addi¬ 
tional meta ls (as in a C C SN). 

From |Groh et al.1 d2014l) . a 60 M 0 star (as inferred by 
IFolev et ali 201 ll for the progenitor of SN 2009ip from pre¬ 
explosion images) will have a final mass of 12 M 0 at the point 
of explosion. Such a star will produce a few solar masses of 
O during hydrostatic burning before it undergoes core-collapse 
dWooslev & Hegen 20071), which is on the order of the ejected O 
mass in SN 1987A Li & McCravl 19921) . If we assume that the den¬ 
sity in SN 2009ip with respect to SN 1987A will scale according 
to the inverse cube of their relative ejecta velocities at late times, 
then the O density in SN 2009ip would be a factor ~ 10 higher than 
the latter. As optical depth evolves with r 2 , one would have to wait 
VTo < ~3) times longer for the transition to optically thin conditions 
in SN 2009ip. As this transition occurs after a few hundred days 
in most core-collapse SNe, it may take several years in SN 2009ip. 
Our latest spectrum of SN 2009ip was taken at ~2 years after its 
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possible core-collapse in August 2012, so the appearance of [O i] 
lines would not be implausible at this epoch. 

The ejected 56 Ni mass from SN 2009ip is another poorly con¬ 
strained parameter. While energy from CSM interaction contributes 
to an unknown fraction of the flux in the late-time lightcurve, it is 
only possible to set an upper limit of the Ni mass. Still, the limit 
derived (M;v, <0.02 M 0 ) seems lower than would be expected for a 
~60 M q progenitor dNomoto et alj 20061) without fallback. 

While [O ij is usually associated with CCSNe, there is a prece¬ 
dent for a stellar object dis playing these lines, in the case of th e 
B supergiant CPD-52°9243 dSwingj[T98il:IWink1er & Wolfll 1 9891) . 
Interestingly, CPD-52°9243 is classified as a B[e] star; a class of 
massive stars which show evidence for a circumstellar disk, rem¬ 
iniscent o f the torus around SN 2009ip inferred from spectropo- 
larimetry (Mauerhan e t alj|2014h . We also note that the appearance 
of [O i] 4/16300,6364 is not ubiquitous in purported SNe - for exam¬ 
ple, no strong [O i] was seen in the Type Iln SN 1995G at ~2 years 
from explosion {Pastorello et alj2002j ). while in contrast these lines 
were clearly seen for SN 1995N {Fransson et all2Qo3) . The reason 
for the absence of the lines in some SNe is unclear; in some cases 
it may indicate that these are not genuine CCSNe but are rather 
SN impostors, while in other cases (such as the bright, and slowly 
declining SN 1995G) it may simply reflect higher densities. 


SN lines is evidence against core-collapse, however the density in 
the ejecta may still be too high for their formation. 

Finally, the physical mechanism powering the outbursts of 
SN 2009ip prior t o 2012 remains unkn own. As di s cussed in 
Fraser et al. ] l2013ch for SN 201 lht, and in lSmith~et Id] {ml) for 
SN 2009ip, the timing of the late time eruptions so soon before 
a possible core-collapse may suggest a link between the two. The 
availability of well-sampled pre-discovery light curves for a large 
sample of S Ne can help shed light on the ubiquity of these erup¬ 
tions dOfek et ah]|2014l) . and help distinguish between som e of the 
propo sed explanations (e.g|puataert & Shiode 2012; Woosle v~et al] 
120071) for the mass loss. Of course, if SN 2009ip can be conclu¬ 
sively demonstrated to have not undergone core-collapse (for ex¬ 
ample, if it is seen to undergo a new outburst), then the classifica¬ 
tion and interpretation of observationally similar transients such as 
SN 2010mc dOfek et alj|2013al) will necessarily need to be revis¬ 
ited. 

Continuing observations of SN 2009ip are required to moni¬ 
tor the photometric and spectroscopic evolution of SN 2009ip; and 
in particular to check if it fades significantly below the magnitude 
of the progenitor candidate. Alongside this, detailed spectral and 
hydrodynamic modelling is required to interpret the wealth of ob¬ 
servational data obtained so far, to securely identify the lines seen 
in SN 2009ip, and to form a clearer picture of the true nature of SN 
2009ip. 


6 CONCLUSIONS 


Questions remain as to if, how, and why SN 2009ip exploded. 
There is now evidence that some very massive stars can explode 


as Ty pe Iln SNe during a luminous blue variable (LBV) phase 
(e.g. lTrundle et all2008 : lGal-Yam & Leonardl2009l) . According to 


stellar evolutionary theory, such massive stars are not expected to 
explode as H-rich LB Vs, but rather are believed to lose their en¬ 
velo pes to become Wolf-Rayet stars, before exploding as Type Ibc 
SNe. lGroh et all ( 120131 . 120141) have recently used models of massive 
stars, coupled with a model atmosphere code, to predict the appear¬ 
ance of massive stars shortly prior to collapse. From detailed mod¬ 


elling of a rotating 60 M 0 star (similar to the progenitor proposed 
for SN 2009ip bv IFolev et alj|201 il) . they find that it will undergo 
an LBV phase between the end of core H-burning and the start of 
core He-burning. The star will ultimately evolve to a WO star, and 
explode a H-poor SN. While the peripheral location of SN 2009ip 
within NGC 7259 is suggestive of a low metallicity progenitor, it is 
unclear whether this could reduce the mass loss rate by enough for 
a 60 M 0 star to retain its H envelope to the point of core-collapse. 
[Heger et~aD {2003h find that while their non-rotating models above 
~40 M q result in a H-free SN, at lower metallicities the star will 
collapse directly to form a black hole, with no accompanying opti¬ 
cal display. 

Neither the core-collapse nor the non-core-collapse scenarios 
can be entirely excluded yet for SN 2009ip. The energetics and 
ejecta velocities seen during the 2012a and 2012b events are closer 
to those of a core-collapse SN than a SN impostor. On the other 
hand, broad (albeit_much_weaker) lines were seen in SN 2009ip as 
early as 2011 {Pastorello et al.ll2013l) . Furthermore, the low ejected 
Ni mass is difficult to reconcile with the collapse of a ~60 M e star, 
unless fallback is invoked. However, even fallback presents prob¬ 
lems, as one would presumably have to fine-tune a model to hide 
nearly all the Ni and core ejecta from SN 2009ip, but still obtain 
ejecta velocities of ~ 8500 kms -1 as seen during the 2012a event. 
It could also be argued that the apparent absence of broad nebular 
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Figure 18. The computed emission flux from a model of the inner H zoneof a c ore-c ollapse SN, showing only line flux arising from material present in the 
envelope of the pre-explosion progenitor, based on the model of ljerkstrand et all 12012 ) at +400 days. Also shown is the spectrum of SN 2009ip obtained on 
2014 October 2, approximately 400 days after the start of the 2012a event. 
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